Search for VHE gamma rays from SS433/W50 with the CANGAROO-II telescope 
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, ^ , SS433, located at the center of the supernova remnant W50, is a close proximity binary system consisting of a compact star and 
a normal star Jets of material are directed outwards from the vicinity of the compact star symmetrically to the east and west. 
Non-thermal hard X-ray emission is detected from lobes lying on both sides. Shock accelerated electrons are expected to generate 
VHE gamma rays through the inverse-Compton process in the lobes. Observations of the western X-ray lobe region of SS433/W50 
system have been performed to detect VHE gamma rays using the 10m C ANGAROO-II telescope in August and September, 200 1 , 
and July and September, 2002. The total observation times are 85.2 h for ON source, and 80.8 h for OFF source data. No significant 
excess of VHE gamma rays has been found at three regions of the western X-ray lobe of SS433/W50 system. We have derived 99% 
confidence level upper limits to the fluxes of gamma rays and have set constraints on the strengths of the magnetic fields assuming 
the synchrotron/inverse-Compton model for the wide energy range of photon spectrum from radio to Te V. The derived lower limits 
^\ are 4.3 juG for the center of the brightest X-ray emission region and 6.3 ftG for the far end from SS433 in the western X-ray lobe. In 
addition, we suggest that the spot-like X-ray emission may provide a major contribution to the hardest X-ray spectrum in the lobe. 

Key words: gamma rays: observations, ISM: individual (W50), jets, stars: individual (SS433) 



> 

m 
m 



o 



^ PACS: 95.85.Pw, 98.38.-j 



1. Introduction 

The galactic SNR W50 is a strong non-thermal radio source. 
Radio images of W50 show a structure extended over ~2° x 
1 ° with limb-brightened "ears" at the eastern and western ends 
Eiaii. The radio emission from W50 peaks at 419 mJy at 
4.75 GHz, with the spectral index varying from a = 0.3 to 1.0 
(where S oc v^") over the source in the frequency range 0.41- 
4.75^GHz fit]. The distance to W50 is estimated to be 5.5 kpc 
fs, 6], and the age is assumed to be about 10"^ years. SS433, lo- 
cated at R.A. (J2000) = 19''ll"'49^ Dec. (J2000) = -h04°58'48" 
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is the jet source located at the center of W50 with a V-band op- 
tical magnitude of 14.2 [f]. SS433 is a close proximity binary 
system with an orbital period 13.1 days, consisting of a compact 
star and a normal star Jets of material are directed symmetri- 
cally outwards from the vicinity of the compact star to east and 
west at a speed of about 0.26c [8J. The axis of the jets is pre- 
cessing in a cone with half-angle of 20°. The precession period 
is ~ 163 days, and the system is oriented at an angle of ~79° 
to the line of sight |8]. The compact star is not yet identified 
whether it is a black hole or neutron star An evaluation is sug- 
gested by Hillwig et al. (2004) [9] that the system consists of 
a low-mass black hole with a mass of 2.9+0.7 M© and a type 
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A3-7 I supergiant with a mass of 10.9+3.1 Mq. 

Bright diffuse X-ray lobes on the eastern and western sides 
of SS433 were discovered by the Einstein observatory in 1983 
ifioll . and were confirmed by ROSAT and ASCA measurements 
in 1997 01 in . These X-ray lobes are believed to have been 
formed by the jets emitted from SS433, as they are symmet- 
rical about SS433 and lie along the same axis as that defined by 
the east-west elongation of the W50 radio shell and are within 
the precession cone of the SS433 jets. Initial ROSAT and ASCA 
measurements of the spectra of these X-ray lobes were com- 
patible with both power-law model and thermal models [11]. 
Subsequently, however, the ASCA team reported harder power- 
low spectra of the western X-ray lobe in 2000. The spectra 
were found to become softer with distance from SS433 within 
the range of the photon index from F = 1.38 to 2.39 il2lll3ll . 



previously. The region "/?2" includes the center of the brightest 
X-ray region and "/?3" is the edge of the brightest region. The 
latter two regions have harder X-ray spectra than "w2". 



They also reported that the thermal model is acceptable only 
when an unusually low metal abundance is assumed. This indi- 
cated that the X-ray emissions of the western X-ray lobe have a 
non-thermal origin. These results were consistent with the pic- 
ture that the high energy electrons are generated at the sides of 
the X-ray lobes closest to SS433, with the electrons losing their 
energy by the synchrotron emission as they travel further from 
SS433. The ASCA results suggest that very high energy elec- 
trons with energies up to several hundred TeV are expected to 
be produced through shock acceleration in the X-ray lobes, and 
the VHE gamma ray emissions are expected to be generated 
through the inverse-Compton (IC) scattering. 

In 2005, the HEGRA team gave the flux upper limits on the 
VHE gamma rays at a few percent of the Crab nebula flux for 
the regions reported by the ROSAT and ASCA team in 1997 |fl4|. 
In an effort to detect VHE gamma rays, we observed the west- 
ern X-ray lobe region with the CANGAROO-11 air Cherenkov 
imaging telescope in 2001 and 2002, based on the ASCA result 
of 2000. We report the results of our observations and discuss 
the possible environmental conditions of the western X-ray lobe 
region. 

2. Observations 

The 10 m CANGAROO-11 telescope f]3\ is located near 
Woomera, South AustraHa (136°47'E, 31°06'S, 160m a.s.l.) 
and consists of 1 14 segmented spherical mirrors each of 80 cm 
diameter [16]. An imaging camera consisting of 552 PMTs 
is placed at the focal plane covering a field-of-view (FOV) of 
2?76x2?76. The CANGAROO-11 telescope has an angular res- 
olution of 0?30 (29 pc at a distance of 5.5 kpc) with an energy 
threshold of 850 Ge V for a Crab-like energy spectrum. 

The ASCA team reported the results of three regions in the 
western X-ray lobe in 2000, based on their high spatial and 
spectral resolving power over a wide energy range il2l 11311 . 



These regions were named positions 1, 2 and 3 by the ASCA 
team in order of distance from SS433, centered 23', 31' and 39' 
west of SS433, respectively. Hereafter, we call these regions 
"pi", "p2" and "p3", respectively (Fig.[T]). According to the 
ASCA results, the region "pi" shows the hardest X-ray spec- 
trum of the three regions and has a harder X-ray spectrum than 
the region "wl" which the ROSAT and ASCA team reported on 
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Figure 1 : 4.75 GHz radio contour map of W50 overlaid on X-ray images of 
tlie western X-ray lobe (dots) obtained by ASCA 1 12]. Solid white circles show 
the regions reported by the ASCA team in 2000. In this paper, we label these 
regions "pi", "p2" and "p3'\ Solid black circles show the regions repoiled by 
the ROSAT and ASCA team in 1997, which are labeled "wl" and "w2". 

Observations of the western X-ray lobe region were per- 
formed using the CANGAROO-11 telescope in August and Septem- 
ber 2001, and in July and September 2002. Based on the ASCA 
results, the tracking position was set at "pi", since the region 
showed the hardest power-law spectrum with the photon in- 
dex of F = 1.38, suggesting shock acceleration could be tak- 
ing place. The data were obtained by ON source observations 
tracking the source position, and by OFF source observations 
for the background. ON source observations were timed to con- 
tain the meridian passage of the target, as was done by Enomoto 
et al. [17]. Thus, the maximum elevation angle during the ob- 
servations was about 54 degrees. OFF source observations were 
timed to have suitable off'set right ascensions which varied day 
by day from -l''36'"48* to -h4''48'"18'. The total observation 
times in 2001 were 5 1 .4 h (ON source) and 49.5 h (OFF source) 
and, in 2002, 33.8 h (ON source) and 31.3 h (OFF source). 

3. Analysis 

We performed a preselection for the data analysis to obtain 
the data under good weather conditions. To reduce night-sky 
background noise, we selected air shower events having at least 
4 adjacent pixels with more than ~3.3photoelectrons, which 
resulted in a stable shower rate. Periods of data with a shower 
rate less than 1 .5 Hz were not used for the present analysis to 
eliminate the effects of partial cloud and dew formation on the 
mirrors. Moreover, in order to stabilize the shower rate, data 
taken at elevation angles less than 40° were not used. After 
these selections, 60.6% of the total observation time was used 
for the analysis (Table [T]i. The resulting mean elevation angle 
was approximately 51°. 

To reduce the effects of the night-sky background, we have 
used the timing information. The pixels which were triggered 
more than 30 nsec from the average trigger time of a shower 
were eliminated. For each pixel, trigger counts within a 700 //i 
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Table 1: Observation time, ta\,^ hours, and selected time, ^s^i hours, in 2001 and 
2002. 
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2001 


36.3 
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51.4 


30.4 


/ 


49.5 


2002 


17.2 


/ 


33.8 


16.7 


/ 


31.3 


Total 


53.5 


/ 


85.2 


47.1 


/ 


80.8 



to increase the accuracy of the orientation angle of the im- 
age, Alpha 123, 21]. To differentiate gamma ray like events 
from cosmic ray like events, we adopted the Likelihood method 
1 2211 which has a higher selection efficiency for gamma rays 
than a conventional parameterization technique. Figure|2]shows 
the MC results for the distributions of Width and Length of 
gamma rays, assuming a spectral index of -2.5. The observed 
points for cosmic rays (background) are shown by the dots in 
the same figure. The Likelihood method uses a single param- 



period were recorded once a second, and were checked during 
the off-line analysis to exclude pixels having high trigger rates, 
which were generally caused by the passage of a bright star 
through the FOV. Trigger counts were summed run by run to 
search for the effect of stars passing through the FOV. After 
excluding pixels having more the 15 triggers per 700 fis period, 
there were no apparent effects of stars during the whole period 
of observations, including the brightest star with the magnitude 
of 4.9 in the OFF source observations on September 6 and 7 
in 2002. Further, in each set of the data, we also eliminated a 
small number of pixels which showed deformed ADC spectra. 
The deformed ADC spectra were determined by the following 
procedure. 

1 . The ADC spectrum of each pixel was made for both years 
using the data which satisfied the trigger condition. 

2. The averaged ADC spectrum of 48 reference pixels which 
were located at symmetrical positions in the focal plane 
with respect to the pixel being examined was defined as 
a reference for each pixel. 

3. ;t'^ of the ADC spectrum against the corresponding refer- 
ence ADC spectrum was obtained c^sx\DC ^'^^ sach pixel. 
In addition, the number of events which satisfied the trig- 
ger condition was compared with averaged number for 
the reference pixels, and the of this number was ob- 
tained as A'entry ^'^ch pixcl, tOO. 

4. Pixels which had larger values ofx^^c than se- 
lected threshold values were eliminated in order to obtain 
good shower images, since these were possibly due to a 
hardware fault. 

After performing this procedure, the pixels having the deformed 
ADC spectra or high trigger rates were eliminated from the 
data of both years. For each shower event, the lower energy 
events which have smaller size shower images tend to be de- 
formed by the hardware noise. To avoid this effect, we selected 
the events above the threshold SUMADC value (sum of ADC 
values of triggered pixels) of 2100 which corresponds to about 
23 photoelectrons. 

The analysis of the data was performed based on the imag- 

El. We calculated 



ing atmospheric Cherenkov technique 11181 
the imaging parameters (Hillas parameters). Distance, Length, 
Width using Monte Carlo (MC) simulations for gamma rays and 
OFF source data for cosmic rays. We selected the events under 
the conditions: the distance of the centroid of the image from 
the center of the FOV was less than 1.3° to eliminate the edge 
effect of the FOV of the camera, and 0.2° < Distance < 1.2°, 





Figure 2: Distributions of Width (left panel) and Length (right panel) for gamma 
ray MC data (contours) and OFF-source experimental data (dots) in the energy 
range 0.2-50 TeV. Each parameter depends on SUMADC which is the sum of 
ADC values event by event, and is approximately proportional to the energy of 
the incident particle. 

eter: Lratio = Py/ {Py + Pcr), where Py and Pcr are the prob- 
ability of the event being due to a gamma ray and a cosmic 
ray, respectively. Both probabilities can be estimated from the 
products of individual probabilities for Width and Length which 
are derived from the probability density functions (PDFs), in- 
cluding its energy dependence. The PDFs were obtained using 
the MC for gamma ray initiated showers and the OFF source 
data for cosmic rays. Figure [3] shows the distributions of Lratio 
expected for gamma rays and cosmic rays. In the region with 
Lmio < 0.35, cosmic rays exceed gamma rays, whereas in the 
region with L,atio ^ 0.35, gamma rays exceed cosmic rays. In 
this analysis, we used the data with Lratio > 0.4 to select the 
candidates of gamma rays. The subsequent selection of events 
with Alpha < 20° eliminates ~9Q% of the cosmic ray events 
but retains ~60% of the gamma ray events. 

To check the feasibility of our observations and the analysis 
procedure, we analyzed Crab nebula data obtained in Decem- 
ber 2000 using the same analysis code. Total analyzed times for 
ON source and OFF source data were 14.9 h and 13.8 h, respec- 
tively. The maximum elevation angle was about 37 degrees, and 
the energy threshold was estimated to be ~2 TeV. The measured 
gamma ray fluxes within the energy range of 2 - 10 TeV are 
shown in Table |2] and Fig. H] Figure |4] also shows power-law 
spectra obtained by H.E.S.S. and MAGIC [ZHH. Using the 
measured fluxes, a fit for a power-law spectrum gives a differen 
tial flux normalization at 1 TeV of (2.9+2.8stat)xlO 
and a power-law index of 2.57+0.59stat- Although the obtained 
spectrum has relatively large statistical errors, the differential 
flux of CANGAROO-II at 4.6 TeV with a power-law spectrum 
showed good agreement with that of H.E.S.S., within 9% and 
MAGIC, within 14% QIH. 



cm ^ sec ' TeV ' 



3 



4. Results 




0,1 0.2 0,3 0.4 0.5 0,6 0.7 0.8 0,9 1 



Likelihood ratio 



Figure 3: Distributions of Likelihood ratio (L,atio) for tlie gamma ray MC data 
(blank area) and OFF-source experimental data (hatched area). The number of 
the MC events are normalized to the OFF-source data. In this paper, we used 
the data with L„tio > 0.4. 



Table 2: Differential fluxes of the Crab nebula. Only statistical errors are 
estimated in this analysis. 



Mean energy 


Differential flux 


(TeV) 


(photons cm"^ sec"') 


2.2+0.2 


(3.7+3.0)xlO-'^ 


2.9+0.3 


(1.8+1. l)xl0-'2 


4.6+0.4 


(6.3+2.9)xlO-'3 


8.6+0.5 


(l.l+0.6)xl0"'^ 



The left panel of Fig. |5] shows the energy dependence of 
the effective area for gamma rays of the region "pi" obtained 
by the MC simulation. The effective area is almost constant at 
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Figure 5: The effective area for gamma rays of the region "pi" is shown in left 
panel. The right panel shows the effective area multiplied by E^--^ . 

~ 1.5x10^ cm- above2TeV. The right panel of Fig. |5]shows the 
energy dependence of the detection efficiency for gamma rays 
of "pi". Thus the energy threshold for gamma rays of "pi" 
is estimated to be 850 GeV from this distribution. The effective 
areas and the energy thresholds for "p2" and "p3" are estimated 
to be the same as for "pi". 

Figure |6] shows the distributions of the image orientation 
angle. Alpha, at "pi", "p2" and "p3" for the combined data 
of 2001 and 2002. The number of OFF source events were 
normalized to the ON source data in the range of Alpha > 30 
degrees. The normalization factor was ~1.1 for all regions. The 
number of excess events was obtained by subtracting the num- 
ber of the OFF source events from the ON source events in the 
range of Alpha < 20 degrees. 
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Figure 4: Differential fluxes of Crab and spectra from H.E.S.S. and MAGIC. 
Only statistical etTors are shown. 
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Figure 6: Alpha distributions of the combined data for "pi" (left), "p2" (mid- 
dle) and "p3" (right). The points with statistical error bars show ON source 
data and hatched histograms show OFF source data. 

The statistical significances of the excesses at "pi", "p2" 
and "p3" were 0.39, -0.11 and -1.0 cr, respectively. Since we 
found no significant excess from the data, we derived the 99% 
confidence level upper limit fluxes, using Helene's method [25J, 
to be 1.5x10"'^ 1.3x10"'^ and 7.9x10"'^ cm"^ gg^-i, respec- 
tively, for VHE gamma rays with E > 850 GeV. These results 
are summarized in Table |3] 
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Table 3: Results of a search for VHE gamma rays with the CANGAROO-II 
telescope from the western X-ray lobe observed by ASCA. The 99% C.L. upper 
limit fluxes are given above the energy threshold (£th = 850 GeV). 



Source 


R.A. 


Decl. 




"s 


■^99% 






+4°57'46" 


39 


0.39 


1.5 


p2 


19h09™44^ 


+4°58'48" 


-12 


-0.11 


1.3 


p3 


19*'09™12'^ 


+4°59'13" 


-97 


-1.0 


0.79 



' Number of Excess events included in Alpha < 20°. 
' Statistical significance in units of standard deviation cr [2(H . 
99% C.L. upper limit flux for E>Eih in unit of 10"'" phcm"" s"'. 



5. Discussion 

According to the ASCA measurements, the X-ray lobes are 
considered to be generated by non-thermal emissions lll2lll3[l . 
The CANGAROO-II observations were carried out based on 
this ASCA result, aiming to detect the VHE gamma rays, but no 
evidence was found for gamma ray emissions above 850 GeV. 
We obtained upper limit fluxes which are useful to constrain the 
parameters of the western X-ray lobe based on the non-thermal 
model through a synchrotron/inverse-Compton (IC) emission 
process. Therefore, we try to understand the wide energy range 
of photon spectrum from radio to TeV using the synchrotron/IC 
model. 

we extracted five data 



From ASCA X-ray data UJ, |13[1 , 
points in the energy region from 0.7 to lOkeV by fitting to 
the power-law spectrum. Similarly, we used the radio data of 
Downes et al. at 1.7, 2.7 and 4.75 GHz fl 0]. The radio and 
X-ray fluxes of "pV, "p2" and "p3" were calculated from the 
intensities of respective energy ranges correcting the FOV to 
the 8' of ASCA. 

As the seed photons for the IC process, we examined the 
possibility of IR radiation in addition to the cosmic microwave 
background (CMB). First, we checked the IR photon field. An 
upper limit to the IR flux is given by Band (1987) Il27ll for the 
optical filament region observed across "/:i2". Using the upper 
limit flux, we have obtained the upper limits of respective en- 
ergy densities as ~0.033 eV cm"'' at 12 /urn and ~0.047eVcm"-' 
at lOOjum. Since we have no information on IR flux at "pi" 
and "p3", we assume the same upper limits for them. The up- 
per limit values for the interval wavelengths are interpolated 
assuming a power-law spectrum. 

Next, we checked optical photon field. Hillwig et al. (2004) 
10] claimed that the normal star in the SS433 system was a 
type A3-7 I supergiant with a mass of 10.9 + 3.1 Mq. This im- 
plies a typical temperature and stellar radius of rs~8000 °K and 
7?s~40/?Q [28]. Assuming the normal star to be a blackbody ra- 
diator, we estimated the energy density to be ~0.067eVcm"^ 
at the peak frequency of 5.0x10'"* Hz. This value of energy 
density negligibly contributes to the IC process compared to 
the CMB. Further more, we evaluated the energy density of 
the optical filament, since "/?2" includes the optical filament 
in the western lobe. From Boumis et al. (2007) Il29(l . we ob- 
tained a energy density of ^10 "* eV cm for Ha line, and the 
same energy density levels for other emission lines. Also from 
Mazeh et al. (1983) Q, we obtained -lO^^eVcm"^ in Nii 



6584 A line. Comparing these energy densities to that of CMB 
(~0.26eVcm"^) coupled with the Klein-Nishina eff'ects, the 
optical photon field provides almost negligible contribution to 
the IC process. Therefore, for all regions we do not take into 
account the optical photon field as the seed photons of IC emis- 
sion. We also checked the possibilities that both the X-ray and 
radio photons of the western X-ray lobe can contribute as IC 
seed photons, and found negligibly small contributions com- 
pared to die CMB. 

For the spectrum of electrons which drive the synchrotron/IC 
process, we assumed £'"''exp (-/i/Zijnax) where y is the power- 
law index, and fimax is the exponential cutoff energy of elec- 
trons. Recent measurements of Hi absorption/emission spec- 
tra and '^CO spectrum toward SS433 support 5.5+0.2 kpc for 
the distance towards the SS433/W50 system, and constrain the 
age of W50 to be younger than 10^ yr [IsT]. We obtained the 
expected fluxes for the synchrotron/IC model adopting a dis- 
tance of D = 5.5 kpc. The ratio between the size of synchrotron 
emission region and that of IC was assumed to be unity. We 
obtained moderately good fits for "p2" and "/?3" on the inter- 
pretation for the wide energy range of the photon spectrum by 
freeing all four parameters. The fit for "/?2" gives the lower 
limit of the strength on the magnetic field, Bmin, to be 4.3 fiG 
with;^f^ - 7.0/5, and the fit for "p3" gives B,nin - 6.3 pG with 

= 15/5. The photon indices of both regions were estimated 
to be ~1.7. For "pi", on the other hand, it is very difficult to 
understand the wide energy range of photon spectrum with a 
single synchrotron/IC emission model since the X-ray emission 
has a very hard spectrum. Hence, we presented the result un- 
der the assumption of a cutoff energy of electrons at "pi" as 
510TeV given by ~280TeV (Ex/lkeWy^ (B/l fiGT^'^ (Ex; 
the energy of X-rays generated by the synchrotron process), us- 
ing the X-ray energy of lOkeV and the strength on the mag- 
netic field of 3 pG. The obtained spectral energy distributions 
for "pi", "p2" and "p3" are shown in Fig.|71 

The fit for "pi" giving /dof = 420/7 clearly indicates 
that a unified interpretation for the wide wavelengths from ra- 
dio to TeV does not work well at "pi". This result leads to 
an alternative interpretation that the X-ray emissions does not 
share the same radiation mechanism with the radio emissions. 
According to the ASCA result, the radial distribution of the hard 
X-rays (3.0 - 10.0 keV) shows a steep peak of ~3' spread at 
near the center of "pi" [12]. In contrast to such a "hot spot" 
feature, the radial distribution of the soft X-rays (0.7 - 3.0 keV) 
does not show a clear peak in "pi", and is similar to that at 
"p2" and "p3" in both energy bands. Thus, the "pi" emis- 
sion region is likely to be a combination of the spot-like hard 
X-ray emission at the center of "pi" together with diffuse X- 
rays similar to "p2" or "p3". A similar morphological feature 
has been found by the recent X-ray measurement with XMM- 
Newton for the eastern X-ray lobe [32]. Based on this combined 
emission model for "pi", it can be suggested that the diffuse 
emission region shares the radiation mechanism from the radio 
to X-ray energy and the hard X-ray spectrum is mainly domi- 
nated by the spot-like emissions. Therefore, we examined the 
synchrotron/IC emission model without the radio data, since 
the diffuse emissions are preferable to provide minor contribu- 
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Figure 7: Synchrotron/IC model fitting to tlie spectral energy distributions at 
"pi" (to p), "p2" (middle) and "p3" (bottom). The X-ray fluxes are from ASCA 
Il4.ri3tl and the radio fluxes are calculated from the result of the Eft'elsberg 100- 
m telescope for each frequency [ l, 4]. The 99% C.L. upper limit fluxes of this 
work are shown by down-arrows. Their mean energies are 0.9, 1.2 and 3.9 TeV, 
respectively. The thick solid lines represent the synchrotron/IC spectra. Two 
components (CMB and IR) of IC ai'e shown by the dashed hnes and the dotted 
lines, respectively. 



tion for the X-ray spectrum. The alternative spectral energy 
distribution for "pi" is obtained as Fig. [8] assuming the same 
cutoff energy of iSmax - 5 10 TeV and B - 3 fiG as in Fig.|2l 
The resulting parameters for "pi", "p2" and "/?3" are summa- 
rized in Table, m As shown in Fig. [8] the expected flux of IC 
emissions at TeV region does not exceed our upper limit flux, 
hence it is found that the interpretation by this model is accept- 
able. Although this model is a possible case, it seems to re- 
quire the combination of multiple emission components to un- 
derstand the wide energy range of photon spectrum using the 
synchrotron/IC model for "pi". 



[ 1 1 1 1 1 1 1 1 1 1 1 1 

r Synchrotron/IC {E^^_=510TeV) 


1 1 1 1 1 1 ' 


B=1(iG/r=0.96 
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Figure 8: Alternative spectral energy distributions for "pi". The data shown 
here are the same data with Fig. although only the X-ray data are taken into 
account for this fit. The cutoft' energy of electrons is set as same value of £niax = 
510TeV as Fig.0 The solid line, the dashed line and the dotted line represent 
the synchrotron/IC spectra under the various magnetic field assumptions; 6=1, 
3 and 10 /.(G, respectively. 



Region 


^min 

(mG) 


F 

^max 

(TeV) 


"r 




>1 
p2 
p3 


"3.0 
4.3+0.1 
6.3±0.3 


"510 
440+60 
130+10 


1.1+0.0 
1.7±0.0 
1.7±0.0 


3.5/(5-2) 
7.0/(9-4) 
15/(9-4) 



" Lower limit on the strength of the magnetic tield. 

Photon index (=(y + 1) /2). 
^ Result of a fit without the radio data. 
^ Assumed parameters. 



6. Conclusion 



According to the analysis of the ASCA data 0121 11311 for 
SS433/W50 system, the X-ray spectra of 3 regions in the west- 
ern part of the X-ray lobe can be explained by non-thermal 
emission. Therefore, shock acceleration may be present in this 
region and the emission of the VHE gamma rays may result. 

Using the 10 m CANGAROO-II telescope, we have searched 
for gamma rays in the VHE region from the western part of the 
X-ray lobe of SS433/W50 system. We detected no significant 
excess of gamma rays from this region. To check the reliabil- 
ity of our observations and the analysis procedure, we analyzed 
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Crab nebula data using the same analysis code, and obtained 
consistent results with recent measurements by H.E.S.S. and 
MAGIC. The 99% confidence level upper limits on the fluxes 
of gamma rays for '"pi", '"p2" and "/>3" were as 1.5x10"'^, 
1.3x10-12 and 7.9x10-1^ cm-^ sec"! above 850GeV, respec- 
tively. Using these upper limit fluxes, we derived the lower 
hmits of the magnetic field as to be 4.3 and 6.3 yuG for "p2" and 
"/?3", respectively, under the assumption of a synchrotron/inverse- 
Compton model for the wide energy range of photon spectrum 
from radio to TeV. The same interpretation for "pV was at- 
tempted and found to be diflftcult. However, we suggested the 
alternative interpretation for "pV by assuming the combined 
X-ray emissions which consists of the difliised X-ray emissions 
and the spot-Uke hard X-ray emissions. Since the spot-like 
emissions were supposed to provide major contribution to the 
hard X-ray spectrum, we examined the model without the radio 
data and found this interpretation was acceptable. 
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